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We update the ground-state energy ratio of unitary Fermi gas to noninteracting Fermi gas (^) 
from the e expansion by including the next-to-next-to-leading-order (NNLO) term near two spatial 
dimensions. Interpolations of the NNLO e expansions around four and two spatial dimensions with 
the use of Fade approximants give ^ ~ 0.360 ± 0.020 in three dimensions with the uncertainty due 
to different interpolation functions. This value is consistent with the previous interpolations of the 
NLO e expansions ^ ~ 0.377 ± 0.014 in spite of the large NNLO corrections. 
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I. INTRODUCTION 

Two-component fermions interacting via a zero-range 
and infinite scattering length interaction have attracted 
intense attention across many subfields of physics [1] . Ex- 
perimentally, such a system can be realized in trapped 
atoms using the Feshbach resonance and has been exten- 
sively studied [2]. The most important property of the 
system is the scale invariance of the interaction, and thus, 
it can be thought of a rare realization of nonrelativistic 
conformal field theories [3-6] . 

As a consequence of the scale invariance of the interac- 
tion, all physical quantities at finite density and zero tem- 
perature are determined by simple dimensional analysis 
up to dimensionless constants of proportionality. Such di- 
mensionless parameters are universal depending only on 
the dimensionality of space. A representative example 
of the universal parameters is the ground-state energy of 
the Fermi gas at infinite scattering length (unitary Fermi 
gas) normalized by that of a noninteracting Fermi gas 
with the same density: 

(1) 

Here we put a subscript d to emphasize that S^a is a func- 
tion of the dimensionality of space. Because is a fun- 
damental quantity characterizing the unitary Fermi gas, 
there have been substantial efforts to determine its value 
in c? = 3 both from experiments [7-13] and Monte Carlo 
simulations [14-23]. 

For analytical treatments, the scale-invariant interac- 
tion implies great difficulties because there seems to be 
no parameter to control a theory. However, it was shown 
that the problem of unitary Fermi gas can be solved sys- 
tematically with appropriately formulated perturbation 
theories if the dimensionality of space d is close to 4 or 
close to 2 [24-26]. This is inspired by the special na- 
ture of four and two spatial dimensions for the zero-range 
and infinite scattering length interaction [27] : the unitary 
Fermi gas becomes a noninteracting Bose gas in d — A 
(^d^4 — *■ 0) , while it becomes a noninteracting Fermi gas 
in d = 2 (fd^2 !)■ Corrections to S^d near four and 
two spatial dimensions have been computed up to next- 
to-next-to-leading order (NNLO) in terms oi e — A — d 



and e = d - 2 [24, 25, 28, 29]: 

(6-e)/(4-e) 

a-. = 2 

X [1 - 0.04916 e - 0.95961 + 0{€^)\ 

and 

C2+e = l-e + 0.80685 e2-hO(e^). (3) 

Because NNLO corrections turn out to be large, naive 
extrapolations of the e and e expansions to the physical 
case in c? = 3 do not work at all. The more appropriate 
way to obtain the value of in d = 3 is to interpolate 
the two expansions. This procedure has been carried 
out by using the next-to-leading-order (NLO) expansions 
around d = 4 and d = 2 [25] and by using the NNLO 
expansion around d — A and the NLO expansion around 
d = 2 [28], and reasonable agreement with results from 
Monte Carlo simulations was found. 

The main purpose of this paper is to update in d = 3 
by including the NNLO term near two spatial dimensions. 
First, we review the interpolations of the NLO e expan- 
sions to see the stability of the results to the choice of in- 
terpolation schemes (Sec. H). We then show results from 
the interpolations of the NNLO e expansions in Sec. HL 
Finally, a summary and concluding remarks are given in 
Sec. IV. The NNLO correction to near d = 2 shown 
in Eq. (3) is computed in the Appendix. 



II. INTERPOLATIONS OF NLO EXPANSIONS 

In order to see the stability of the results to the choice 
of interpolation schemes, we review the interpolations of 
the NLO e expansions by using Fade approximants with 
and without applying the Borel transformation. 



A. Pade interpolation 

The simplest way to interpolate the two expansions 
around d = 4 and d = 2 is to use the Fade approximants. 



FIG. 1: (Color online) The universal parameter as a function of spatial dimensions d. The upper curve is the extrapolation 
from the NLO expansion around d = 4 in Eq. (2), while the lower line is the extrapolation from the NLO expansion around 
d = 2 in Eq. (3). The middle three curves show the Pade (left panel) and Borel-Pade [25] (right panel) interpolations of the 
two NLO expansions. The symbol at d = 3 indicates the result ^3 ~ 0.40(1) from the latest Monte Carlo simulations [21, 23]. 



We write in Eq. (2) in the following form: 



B. Borel-Pade interpolation 



(6-e)/(4-e) 

^4-. = ^ F{e) 



(4) 



where F(e) is an unknown function having the expan- 
sion F(e) = 1 - 0.04916 e - 0.95961 + C»(g3) [35], -yVe 
approximate F{e) by a ratio of two polynomials (Pade 
approximant), 



F[M/N] (e) - 



Po +Pie H KPA/e 



M 



1 + qie 



N 



(5) 



and determine the unknown coefficients so that has 
the correct expansions around d — A and d = 2. If one 
truncates the e and e expansions at NLO, we have four 
known terms and thus Pade approximants F^m/n] satis- 
fying M + N = 3 are possible. We exclude the possibility 
of F[2/i]{e) because it has a pole in a range < e < 2, 
while we expect a smooth behavior of as a function of 
2 < d < 4. 

The left panel in Fig. 1 shows the universal parameter 
£,d as a function of d. The middle three curves show the 
Pade interpolations of the two NLO expansions with the 
use of -F[3/o], F[i/2], and ^"[0/3]- In d = 3, these interpo- 
lations, respectively, give 



^3 « 0.391, 0.366, 0.373. 



(6) 



These three values have an average 0.377 and span a 
small interval ^3 « 0.377 ±0.014. We note that the same 
interpolation scheme was employed to compute the low- 
est two energy levels of three fermions in a harmonic 
potential, and excellent agreement with the exact results 
was found in arbitrary spatial dimensions 2 < d < 4 [5]. 



The other way to interpolate the two expansions is 
to apply the Borel transformation and then use the Pade 
approximants [25]. We first rewrite the unknown function 
F{e) in Eq. (4) in the form of the Borel transformation: 



1 1'°° 

F(e) = - / dte-*/'G{t). 
e Jo 



(7) 



c„e'% the Borel 

oc 



If F{e) has an expansion F{e) = J2n=o ^n. 
transform G(t) has an expansion G(t) — X^J^o 
and thus, the Borel transformation makes the expan- 
sion faster convergent. Then we approximate G(t) by 
the Pade approximant. 



G 



\MIN\ 



M 



PO +pit^ \-pMt 

1 + gitH + qNt'^ 



(8) 



and determine the unknown coefficients so that has 
the correct expansions around d = 4 and d = 2. 

The right panel in Fig. 1 shows the universal parameter 
as a function of d. The middle three curves show 
the Borel-Pade interpolations of the two NLO expansions 
with the use of G[3/o], ^[1/2], and G[o/3]. The possibility 
of G[2/i] is excluded because we could not find a solution 
satisfying the constraints of Eqs. (2) and (3). In d = 3, 
these interpolations, respectively, give [25] 



~ 0.391, 0.364, 0.378. 



(9) 



These three values have an average 0.378 and span a 
small interval ^3 « 0.378 ±0.013. We note that the result 
of Gp/o] is equivalent to that of ^"[3/0] in Eq. (6). 

Comparing the results in Eqs. (6) and (9), one can 
see that the interpolated values do not depend much on 
the choice of the Pade approximants and also the Borel 



3 




FIG. 2: (Color online) The universal parameter as a 
function of spatial dimensions d. The lower right curve is 
the extrapolation from the NNLO expansion around d — 4 
in Eq. (2), while the upper left curve is the extrapolation 
from the NNLO expansion around d = 2 in Eq. (3). The 
middle four curves show the Fade interpolations of the two 
NNLO expansions. The symbol at d = 3 indicates the result 
^3 « 0.40(1) from the latest Monte Carlo simulations [21, 23]. 



transformation does not improve the interpolated val- 
ues. The latter can be understood because we have only 
a few terms in the expansion over e (two terms up to 
NLO), and thus the advantage to apply the Borel trans- 
formation is little. This situation does not change even 
if we include the NNLO term near c? = 4. The use of 
the Borel transformation may become important once we 
have more higher-order corrections and the e expansion 
is not a convergent series. 

III. FADE INTERPOLATION OF NNLO 
EXPANSIONS 

We now include the NNLO terms near d = 4 and d = 2 
to interpolate the two expansions. Here we only use the 
Pade interpolation by the above-mentioned reason [36]. 
Because we have six known terms, Pade approximants 
F[M/N] satisfying M -I- = 5 are possible. However, we 
exclude the possibility of ^"[2/3] and i^[i/4] because they 
have poles in a range < e < 2, while we expect a smooth 
behavior of £,d as a function of 2 < rf < 4. 

Figure 2 shows the universal parameter S,d as a function 
of d. The middle four curves show the Pade interpola- 
tions of the two NNLO expansions with the use of /q] j 
i^[4/i], ^"[3/2], and ^"[0/5]- In c? = 3, these interpolations, 
respectively, give 

^3 « 0.340, 0.372, 0.370, 0.357. (10) 

These four values have an average 0.360 and span an 
interval ^3 « 0.360 ± 0.020 [37]. It is understandable 



that the interpolations of the NNLO expansions have the 
larger uncertainty because of the large NNLO corrections 
both near d — 4 and d = 2 [see Eqs. (2) and (3) and also 
Fig. 2] . What is remarkable is that in spite of such large 
NNLO corrections, the interpolated values are consistent 
with the previous interpolations of the NLO expansions 
^3 ~ 0.377 ± 0.014. Therefore we conclude that the in- 
terpolated results are stable to inclusion of higher-order 
corrections and thus the e expansion has a certain pre- 
dictive power even though the knowledge on higher-order 
terms in the expansions over e = A — d and e = d — 2 is 
currently lacking. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this paper, we have updated the ground-state energy 
ratio of unitary Fermi gas to noninteracting Fermi gas 
(^) from the e expansion by including the NNLO term 
near two spatial dimensions. We found that the Pade 
interpolations of the NNLO expansions around d — A and 
d = 2 give £, « 0.360 ± 0.020 in d = 3 with the relatively 
small uncertainty from different interpolation functions. 
Although the NNLO corrections are large both near d = 
4 and d ~ 2, the interpolated value is consistent with 
the interpolations of the NLO expansions ^ « 0.377 ± 
0.014. This indicates that the interpolated results are 
stable to inclusion of higher-order corrections, and thus 
the e expansion has a certain predictive power. Indeed, 
our interpolated values reasonably agree with the results 
from the latest Monte Carlo simulations, ^ w 0.40(5) [21] 
and C < 0.40(1) [23]. 

Our analysis also implies that in order to obtain ap- 
propriate results from the e expansion, it is necessary 
to incorporate the expansions both around d = 4 and 
d = 2. Other than ^ studied in this paper, interpolations 
of NLO expansions around d = 4 and d = 2 have been 
employed to estimate the critical temperature Tc [26], 
thermodynamic functions at [26], and the ground- 
state energy of a few fermions in a harmonic potential [5] . 
Quasiparticle spectrum [24, 25], atom-dimer and dimer- 
dimer scatterings in vacuum [30], the phase structure of 
polarized Fermi gas with equal masses [25, 31] and un- 
equal masses [29], BCS-BEC crossover [32], momentum 
distribution and condensate fraction [29] , low-energy dy- 
namics [33] , and energy-density functional [34] have been 
studied only in the expansions over e = 4 — d. It is pos- 
sible to obtain better understanding of these subjects by 
further incorporating the expansions in terms of e = d— 2. 
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APPENDIX: NNLO CORRECTION TO NEAR 

d = 2 

In this appendix, we briefly review the e expansion for 
the unitary Fermi gas around two spatial dimensions and 
compute the NNLO correction to £,d in terms of e = d~2 
shown in Eq. (3). The detailed account of the e expansion 
is found in Ref . [25] . 

1. Lagrangian and power counting rule of e 



--•►--(^^^^-•►-- + = 0(e) 

(a) (b) 

FIG. 3: Power counting rule of e. The self-energy diagram 
of ip field (a) is combined with the vertex from £2 (b) to 
achieve the simple e counting. Solid (dotted) lines represent 
the fermion (auxiliary field) propagators iG (iD). 



The unitary Fermi gas near two spatial dimensions is 
described by the sum of following Lagrangian densities 
(here and below h = 1): 

£0 = E + (A.l) 



£2 = V*^- (A. 3) 

Here we have neglected the condensate (/)o ~ /ie~^/'^, 
because its contribution is negligible compared to any 
power corrections of e. 

The first part £0 generates the propagator of fermionic 
field '(/'crj 



1 



Po 



(A.4) 



where £p = /{2m) is the kinetic energy of nonrelativis- 
tic particles. The second part £1 describes the interac- 
tion between fermions mediated by the auxiliary field ip. 
The first term in £1 gives the propagator of <y9, 



(A.5) 



and the last two terms give vertices coupling two fermions 
with ifi. The coupling constant g is given by 



9 



V"2^/ 



(A.6) 



Here the factor (m/i/27r)~'^^'* was introduced so that the 
product of auxiliary fields ip*ip has the same dimension 
as the Lagrangian density. We emphasize that the choice 
of this factor is arbitrary, if it has the correct dimension, 
and does not affect final results because the difference 
can be absorbed by the redefinition of ip. The particular 
choice of g in Eq. (A.6) will simplify expressions for loop 
integrals in intermediate steps. 

If we did not have the last part £2, we could integrate 
out the auxiliary fields if and 1^9*, which leads to 



(A.7) 



which represents the contact interaction of fermions with 
the small coupling g^ ~ e. Therefore, the unitary Fermi 



gas near two spatial dimensions is simply described by 
a weakly interacting system of fermions. The vertex in 
£2 plays a role of a counterterm so as to avoid double 
counting of a certain type of diagrams which is already 
taken into £1. 

The power counting rule of e is summarized as follows. 

1. For any Green's function, we write down all Feyn- 
man diagrams using the propagator from £0 and 
the vertices from £1. 

2. If there is any subdiagram of the type in Fig. 3(a), 
we add the same Feynman diagram where the 
subdiagram is replaced by the vertex from £2 in 
Fig. 3(b). 

3. The power of e for the given Feynman diagram is 
simply 0(e^s/2)^ where Ng is the number of cou- 
plings g. 

Here the dimensional regularization of loop integrals is 
assumed. 



2. Computation of the pressure 

The pressure of unitary Fermi gas has been computed 
up to the next-to-leading order in e [25]. To the lead- 
ing order, the pressure is given by that of nonintcracting 
fermions: 



Pfvcc 2 



dp 



ill ~ Ep) 0{^1 - £p) 



2^ /TO^\'^/2 



(A., 



2) V 27r y 



The next-to- leading-order correction is 0(e), which cor- 
responds to the mean-field correction 



dp 



1 2 



e/i /m/^\'='/2 



(A.9) 



W^\2nJ 



To the next-to-next-to- leading order in e, the pressure 
receives O(e^) corrections from two three-loop diagrams 
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FIG. 4: Vacuum diagrams contributing to the pressure to 
the next-to-next-to-leading order in e. In the right diagram, 
the counter vertex in Fig. 3(b) for each bubble subdiagram is 
implicitly understood. 



depicted in Fig. 4. The left diagram is easily evaluated 
as 



Psa — 



2 



r(i)r(i + i) ^27ry 

Now the right diagram in Fig. 4 is written as [29] 



(A.IO) 



_2 / dkdp 

Psb — g 



!2n 



1 + 9' 



2d 



dq 



-q+ 



(27r)" 



(A.ll) 



where the frequency integrations are already performed. 
We note that -f 1 in the square brackets comes from the 
counter vertex in £2- Due to the 6 functions, the ranges 
of integrations over Sk, £p, and Sq are limited to < £fc < 
4/1, < £p < Ap, and A^ < Sq, where 



Ap = 



(A.12) 



and 



cosxql y/£k + J 'if^- Sk sin Xq 



with cosXp = k ■ p and cosxg = k ■ q. The integration 
over Eq can be performed analytically using dimensional 
regularization. As a result, the expression in the square 
brackets in Eq. (A.ll) becomes 



[•••] 



2 2 



0(e2). (A.14) 



Then, introducing dimensionless variables z = Sk/fJ-, 



pig) 



Ap(q)//i and performing the integration over 



Ep/ii, we obtain the following expression for P; 



P: 



3b 



at- 

'0 JO ^ Jo ^ 
X jA^lnA, -(A, -Ap)ln(A, -Ap)-Ap| 



(A.15) 

Finally the numerical integrations over z, Xpi ^iid Xg l^ad 
to 

2 

Psb = -e'^ (l + 0.0568528) + ©(e^). (A.16) 

Consequently, we obtain the pressure up to the next- 
to-next-to-leading order in e as 



P = Ph. 



P2 + Paa + Pah 



Pfree [l + £ + 0.6931472 1^ + 0{f)\ 



(A.17) 



The universal parameter of the unitary Fermi gas in 
Eq. (1) can be equivalently expressed as = m/^f- From 
the thermodynamic relationship n = dP/dfi and the def- 
inition of the Fermi energy in d spatial dimensions. 



27r 

ep = — 
m 



2/d 



we can determine S.d from the e expansion to be 

- [1 + f- - 0.6931472 

= l-e-|-0.8068528e2-f 0(e3). 



(A.18) 



(A.19) 



(A. 13) This is the result shown in Eq. (3). 
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